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1. Introduction 

Discrete quantum systems were studied originally by and Schwinger!^ and later 

by many authors.^ ^ However, many concepts that appear sharp for continuous systems 
become fuzzy when one tries to apply them to discrete ones. The reason is that the in 
the continuum we have only one harmonic oscillator, while for finite systems, there are a 
number of candidates for that role, each one with its own virtues and drawbacks.^ 

Coherent states constitute an archetypical example of the situation: for the standard 
harmonic oscillator they are well understood, and sensible generalizations have been de- 
vised to deal with systems with more general dynamical groups.'^ However, in spite of 
interesting advances, " the discrete counterparts are still under heated discussion. 

Number states are eigenstates of the Fourier transform, but they are not the only ones: 
coherent states, with their associated Gaussian wave functions, also are. In our opinion, 
this subtle, yet obvious, observation has not been taken in due consideration in this field. 
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From this perspective, a decisive step was given by Mehta,^ who obtained the eigenstates 
of the finite Fourier transform. The purpose of this paper is to further explore this path, 
showing how we can define physically feasible coherent states for n qubits that are also 
Fourier eigenstates. 

Since the natural arena of these discrete coherent states is the phase space, we use the 
discrete Wigner function to picture the corresponding results. The problem of generaUzing 
the Wigner function to finite systems has a long story. Using the notions presented in the 
comprehensive review TO) we construct a Wigner function for these coherent states and 
discuss some of their properties. 



2. Coherent states for n qubits 

We consider n identical qubits (i.e., n noninteracting spin 1/2 systems). We recall that the 
Dicke states, belonging to the symmetric subspace of the representation of SU(2)®", are 
given by 



I"' = V , E ^fcd^l' • ■ • . Ifc: Ofc+1, . . . , 0„)) , (1) 



'k\{n-k)\ 

k 

where {Pk] denotes the complete set of all the possible permutations of the qubits. These 
states can be expressed in terms of the elements of the Galois field GF(2") using the 
standard decomposition in the self-dual basis (a short review of the concepts of finite fields 
needed in this paper is presented in the appendix) 

111, I2, . . . , Ife, Ofc+i, . . . , 0„) i-^ llcTi + 1(72 + • • • + l^Tfc + Offfe+i + . . . + Oct„) . (2) 
Now, consider a SU(2) coherent statd^ 



1 / ! 

\^)= ^l^\^\^Yl. ]:J^^,e\n.k), (3) 

where the complex number ^ is related with the angular coordinates {d, ip) on the Bloch 
sphere by 

£, = cot(?9/2) e~"^ . (4) 
Using the previous correspondence, |^) can be recast as 

10 = , 1' y./2 E ^'^^^W- (5) 

where the function ^1(7), when applied to the field element 7 = 2fe=i 7fe'^fe indicates the 
number of nonzero coefficients 7^ . 
In this case, the Fourier operator is 

^ = ^ E (6) 

/i,i^GGF(2") 
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X being an additive character defined in ( |A.4| i. As it is well known, = 1, so if we 
impose that the states |^) are also eigenstates of F we are lead to 



m = ±10 



(7) 



This immediately implies (all the spins are pointing in the same direction) that there are 
two SU(2) coherent states (with ^± = ±\/2 — 1) that simultaneously are eigenstates of 
the Fourier operator: they are precisely our candidates to be coherent states for n qubits. In 
particular, satisfy the following condition 



1 



2"/2 



Ai,7eGF(2") 



AieGF(2") 



or, equivalently, 



2"/ 



(8) 



(9) 



7eGF(2" 



and the minus sign may appear only for odd number of qubits. 

Equation (|5]l is the abstract form of the SU(2) coherent state. It factorizes in a product 
of single-qubit states when represented in the self-dual basis, i.e.. 



10 



(1 + |^|2)n/2 



ci,...,c„GZ2 



f.m±mL (10) 



Ck being the expansion coefficients of 7 in that basis. The operator transforming from the 
arbitrary basis {ei, £2, into a factorized form is always a permutation given by 



E (K>...K»((Ml|-..(Mn|) 

AieGF(2") 



(11) 



where 



1=1 1=1 



(12) 



Let us examine the simple yet illustrative example of a two-qubit coherent state. In its 
abstract form it reads as 



10 



1 



1 + 



■(|0>+0'^> + 00+OIO)- 



In the self-dual basis (cr, a^) we have the representation 



|0) = 100) = 






VI 



k) = |io) = 







= 101) = 



1 


Vo/ 



W') = 111) = 



(13) 





' 

Vo/ 

(14) 
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in such a way that 



10 = Y 



^ V 1 



In a non self-dual basis (a, a^) we have 



|0) = |00) = 



and 








, |a) = |10) 







10 = r 



= |oi) 



1 





\a') = 111) 



(15) 





' 
(16) 



(17) 



which cannot be factorized. The transition operator for this case is 

/O 1 0\ 

10 
10 

Vooo i/ 

and it is nothing but a CNOT gate performing the operation 

|00) + |01)^|00) + |11). 



(18) 



(19) 



3. Discrete Wigner function 

To gain further insights into the coherent states |^±) we proceed to picture them in phase 
space. To this end, we first note that, while in the continuous case it is possible to translate 
a state by an infinite distance, this is clearly not possible if the space is finite. To "prevent" 
a state from "escaping" the finite phase space it is natural and convenient to use the field 
GF(2") in the representation of the states. 

In consequence, we denote by |a), with a G GF(2"), an orthonormal basis in the 
Hilbert space of the system. Operationally, the elements of the basis can be labeled by 
powers of a primitive element, and the basis reads 

{|0), |a2"-i = l)}. (20) 

These vectors are eigenvectors of the operators belonging to the generahzed Pauli 
group, whose generators are now defined as 

Z0= Yl X{ap)\a){a\, Xp = ^ \a + l3){a\, (21) 

aeGF(2") a6GF(2") 
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SO that 

Z^Xp = xi.c^P)X0Z^. (22) 

The operators ( |2TI ) can be factorized into tensor products of powers of single-particle 
PauU operators and a^, whose expression in the standard basis of the two-dimensional 
Hilbert space is 

a, = |l)(l|-|0)(0|, (T,-|0)(1| + |1)(0|. (23) 

This factorization can be carried out by mapping each element of GF(2") onto an ordered 
set of natural numbers. As we have akeady seen, a convenient choice for this is the self- 
dual basis, since the finite Fourier transform factorizes then into a product of single-particle 
Fourier operators, which leads to 

Z„ = . . . <T^ , Xf3^Srl^ ®...®al- , (24) 

where (oi , . . . , a„) and (6i , . . . , 6„) are the expansion coefficients of a and (3, respectively, 
in the self-dual basis. 

It was shown that the operators 

D{a,P)=ct){a,P)Z^Xp, (25) 

where (j){a, (3) is a phase, form an operational basis in the discrete phase space. The uni- 
tarity condition imposes the condition (3) = y(— a/3). These displacement operators 
(or phase-point operators in the notation of WootteraUJ) allows us to introduce a Hermitian 
kernel 

^("'/5) = ^ E x{c^^-Py)D(pL,y), (26) 

Ai,i'GGF(2") 

in terms of which we can define a well-behaved Wigner function as 

W^,(a,/3) =Tr[eA(a,/3)], (27) 

where g is the density matrix of the system. 

After some calculations, the Wigner function for our coherent states turns out to be 

^ 1^1 ^_L,v,',eGF{2^) 

(28) 

A plot of this function for the case of three qubits is shown in Fig. 1. We also note that the 
marginal distributions take a very simple form: 

To conclude we wish to mention that it is also possible to introduce the notion of 
1 o 

squeezing for these states . In the basis of the eigenstates of Z^, such an operator has the 
following form 

Sx= E (30) 

KeGF(2") 
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Fig. 1 . Wigner function for a coiierent state of a system of tliree identical qubits. 

The following relations hold 

S\ZaS\ = Z^x-l , S\XaS\ = Xa\ , (31) 

so that 

{i±\slx^Sx\i±) = (e±i^I^aA^5Aie±) • 02) 

4. Conclusions 

In summary, we have formulated a new sensible approach to deal with coherent states for 
a system of n qubits. The associated discrete Wigner function has also been worked out. 
Some related problems, as the behavior under time evolution or the extension to systems 
of qudits, will be addressed elsewhere. 

Appendix A. Galois fields 

We briefly recall the minimum background of finite fields needed to proceed through this 
paper. The reader interested in more mathematical details is referred, e.g., to the excellent 
monograph by Lidl and Niederreiter.tl^! 

A commutative ring is a set R equipped with two binary operations, called addition 
and multiplication, such that it is an Abelian group with respect the addition, and the mul- 
tiplication is associative. Perhaps, the motivating example is the ring of integers Z with the 
standard sum and multiplication. On the other hand, the simplest example of a finite ring 
is the set Z„ of integers modulo n, which has exactly n elements. 

A field F is a commutative ring with division, that is, such that does not equal 1 and 
all elements of F except have a multiplicative inverse (note that and 1 here stand for the 
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identity elements for the addition and multiplication, respectively, which may differ from 
the familiar real numbers and 1). Elements of a field form Abelian groups with respect 
to addition and multiplication (in this latter case, the zero element is excluded). 
The characteristic of a finite field is the smallest integer p such that 

pl = l + l + ... + l = (A.1) 

p times 

and it is always a prime number. Any finite field contains a prime subfield Zp and has 
d = elements, where n is a natural number. Moreover, the finite field containing 
elements is unique and is called the Galois field GF(p"). 

Let us denote as Zp [x] the ring of polynomials with coefficients in Zp. Let P{x) be an 
irreducible polynomial of degree n (i.e., one that cannot be factorized over Zp). Then, the 
quotient space Zp[X]/P(x) provides an adequate representation of GF(p"). Its elements 
can be written as polynomials that are defined modulo the irreducible polynomial P{x). 
The multipUcative group of GF(p") is cyclic and its generator is called a primitive element 
of the field. 

As a simple example of a nonprime field, we consider the polynomial .t^ + x + 1 = 0, 
which is irreducible in Z2. If cr is a root of this polynomial, the elements {0, 1, a, = 
(7 + 1 = cr~^} form the finite field GF(2^) and a is a primitive element. 

A basic map is the trace 

tr(a) =a + a^ + ... + aP , (A.2) 

which satisfies 

tr(a + /3)=tr(a)+tr(/3), (A.3) 
and leaves the prime field invariant. In terms of it we define the additive characters as 



■ tr(Q) 
P 



(A.4) 



Xia) = exp 
and posses two important properties: 

x{a + l3)=x{a)xiP), Yl x(a/3) = p"^o,/3. (A.5) 

aGGF(p") 

Any finite field GF(p") can be also considered as an n-dimensional linear vector space. 
Given a basis {Ok}, (k = 1, . . . , n) in this vector space, any field element can be repre- 
sented as 



a = 

fc=i 



J2^k0k, (A.6) 



with ttfc e Zp. In this way, we map each element of GF(p") onto an ordered set of natural 

numbers a <J4> (oi, . . . , a„). 

Two bases {^i, . . . , 0„} and {9'i , . . . , 6*^} are dual when 

tv{0ke[) = dk,i. (A.7) 
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A basis that is dual to itself is called self-dual. 

There are several natural bases in GF(p"). One is the polynomial basis, defined as 

{l,c7,<72,...,a"-i}, (A.8) 

where cr is a primitive element. An alternative is the normal basis, constituted of 

{<7,<7P,...,C7P""'}. (A.9) 

The choice of the appropriate basis depends on the specific problem at hand. For example, 
in GF(2^) the elements {u, a^} are both roots of the irreducible polynomial. The polyno- 
mial basis is {1, a} and its dual is {cr^, 1}, while the normal basis {a, cr^} is self-dual. 
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